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GEOMETRICAL PHASE OPTICAL ELEMENTS WITH SPACE-VA RIANT 
SUBWAVELENGTH GRATINGS 



FIELD OF THE INVENTION 
5 The present invention relates to optical elements. More particularly it relates to 

geometrical phase optical elements with space-variant subwavelength gratings and 
their uses thereof. 

BACKGROUND OF THE INVENTION 

1 0 The Fancharatnam-Berry phase is a geometrical phase associated with the 

polarization of light. When the polarization of a beam traverses a closed loop on 
the Poincare sphere, the final state differs from the initial state by a phase factor 
equal to half the area Q, encompassed by the loop on the sphere (see S. 
Pancharatnam, Proc. Ind. Acad. Sci. A 44, 247 (1956), and M. V. Berry, Proc. 

15 Roy. Soc. (London) A 392, 45 (1984)). In a typical experiment, the polarization of 
a uniformly polarized beam is altered by a series of space-invariant (transversely 
homogenous) waveplates and polarizers, and the phase, which evolves in the time- 
domain is measured by means of interference (see R. Simon, H. J. Kimble and E. 
C. G. Sudharshan, Phys. Rev. Lett. 61 ,19 (1988), and P. G. Kwiat and R. Y. 

20 * Chiao, Phys. Rev. Lett. 66, 588 (1991)). 

Recently, a Pancharatnam-Beny phase in the space-domain was 
considered. Using space- variant (transversely inhomogeneous) metal stripe 
subwavelength gratings, we demonstrated conversion of circular polarization into 
radial polarization (see Z. Bomzon, V. Kleiner and E. Hasman, Appl. Phys. Lett. 

25 79, 1587 (2001)), and showed that the conversion was accompanied by a space- 
variant phase modification of geometrical origin that effected the propagation of 
the beams (Z. Bomzon, V. Kleiner, and E. Hasman, Opt, Lett. 26, 1424 (2001)). 

Previously Bhandari suggested using a discontinuous spatially varying 
waveplate as a lens based on similar geometrical phase effects (R, Bhandari, Phys. 

30 Rep. 281, 1 (1997)). Recent studies have investigated periodic polarization 

gratings (F. Gori, Opt. Lett. 24, 584 (1999), C. R. Fem^ndez-Pousa, I. Moreno, J. 
A. Davis and J. Adachi, Opt. Lett. 26, 1651 (2001) and J. Tervo, and J. Turunen, 
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Opt. Lett. 25, 785 (2000)). These authors showed that the polarization of the 
diffracted orders could differ from the polarization of the incident beam. We 
intend to prove and utilize a connection between the properties of such 
polarization gratings and the space-domain Pancharatnam-Berry phase. 

5 

BRIEF DESCRIPTION OF THE INVENTION 

There is thus provided, in accordance with some preferred embodiments of the 
present invention a space variant polarization optical element for spatially 
manipulating polarization-dependent geometrical phases of an incident light beam, the 

1 0 element comprising a substrate comprising a plurality of zones of gratings with a 

continuously varying orientation, the orientation denoted by ^(x,^) which is equal to 
half of a desired geometrical phase (DGP) modulus 2tt, each grating with a local 
period that is smaller than the wavelength of the incident light beam. 

Furthermore, in accordance with some preferred embodiments of the present 

1 5 invention, the orientation of the grating satisfies the equation 20ix,y) = ^ V^l^^j^ » 

where/is a desired focal length and A is the wavelength of the incident light beam, 
whereby the optical element is used as a converging lens if the incident light beam 
exhibits right-hand circular polarization, and as a diverging lens if the incident light 
beam exhibits left-hand circular polarization. 

20 Furthermore, in accordance with some preferred embodiments of the present 

invention, the following relation is maintained V x = 0 , where 

= A:o(x,y)[cos(«5>/jc,;/)/2)x + sin(?7rf(x,j/)/2)y], where Kg is a grating vector, 
X and y are unit vectors in the x and y direction, Kq = 27r/A(x^y) , where Ko is the 
spatial frequency of the grating, A is the local period of the grating and (x,y)/2 is 

25 the space-variant direction of the vector so that it is perpendicular to the grating 
stripes at any given point. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the optical element is used as diffraction grating element, wherein the 
following equation is maintained = (l^/d):^^^^ , where d is the period of the 

30 plurality of zones, and wherein the grating satisfies the relation 



wo 2004/003596 



PCT/IL2003/000541 



3 

(^5 y) = / Aq ) sin(;D: / d) exp(-<^ / rf) , where Aq is the local period of the 
grating 2Xy—0, 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the substrate is a wafer. 
5 Furthermore, in accordance with some preferred embodiments of the present 

invention, the wafer is manufactured using photolithography techniques. 

Furfliermore, in accordance with some preferred embodiments of the present 
invention, the wafer is manufactured using etching. 

Furthermore, in accordance with some preferred embodiments of the present 
10 invention, the grating is in blazed form, with opposite blazed directions for incident 
left-hand circular polarization and for right-hand circular polarization, of the incident 
light beam. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the orientation of the gratings varies linearly in a predetermined direction. 
1 5 Furthermore, in accordance with some preferred embodiments of the present 

invention, the orientation of the grating of the zones satisfies the equation 
0{x) = -nxldt^^^ , where ^(x) is the orientation of a grating line at a position x on an 

axis, and d is the period of the orientation of the grating of the zones. 

Furthermore, in accordance with some preferred embodiments of the present 
20 invention, the optical element is used as an optical switch. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the optical element is used as a beam-splitter. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the optical element is used as a Lee-type binary subwavelength structure 
25 mask. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the optical element is used for polarimetry. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the grating in each zone comprise of at least two regions of gratings 
30 arranged in different orientations. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the following relation is satisfied, 0{x^y) = (o{x^y) + c , where x and y are 
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coordinates of a specific position in an orthogonal set of axes, m = arctan(>'/x) is the 
azimuthal angle, and c is a constant. 

Furtheimore, in accordance with some preferred embodiments of the present 
invention, the orientation of the grating is spiral. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the orientation of the grating satisfies the relation = I colly where / is 

a topological charge, and r^co indicate a specific angular position at radius r and angle 

CO, 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the grating satisfy the relation 

o)) = (iTV ro / AoXro / rf""'' cos[(/ / 2 - 1) ] / [/ / 2 - 1] for / ^ 2 , and 
(ry o}) = (2;r ro / Aq ) ln(r / ) for / = 2 , where Ao is the local period of the grating . 

Furthermore, in accordance with some preferred embodiments of the present 
invention, there is provided a space variant polarization optical element for spatially 
manipulating polarization-dependent geometrical phases of an incident light beam, the 
element comprising a substrate comprising a plurality of zones of gratings with 
discretely varying orientation, the orientation denoted by ^(x,>') which is equal to half 
of a desired geometrical phase (DGP) modulus 2Tr, each grating with a local period 
that is substantially smaller than the wavelength of the incident light beam. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the discretely varying orientation comprises rotated orientation. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the rotated orientation varies linearly. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the optical element is used as diffraction grating element, wherein the 
following equation is maintained = (2/r / d)j(\^^^ , where d is the period of the 

plurality of zones, and wherein the grating satisfies the relation 

^gix^y) = (2rf/ Ao)sin(;a:/rf)exp(-^/</) , where Ap is the local period of the 

grating at 3^=0. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the orientation of the grating of the zones satisfies the equation 
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0(x) = -7jx/d\ ^„ , where 0(x) is the orientation of a grating line at a position x on an 

axis, and d is the period of the orientation of the grating of the zones. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the grating in each zone comprise of at least two regions of gratings 
5 arranged in different orientations. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the zones of gratings are arranged in an annular manner. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, the zones of gratings are arranged in a coaxial manner. 
10 Furthermore, in accordance with some preferred embodiments of the present 

invention, there is provided a method for spatially manipulating polarization- 
dependent geometrical phases of an incident light beam, the method comprising: 
providing a substrate comprising a plurality of zones of gratings, with a continuously 
varying orientation, the orientation denoted by ^(jc,;;) which is equal to half of a 
15 desired geometrical phase (DGP) modulus 2tt, each grating having a local period that 
is smaller than the wavelength of the incident light beam 
irradiating the incident light beam onto the substrate. 

Furthermore, in accordance with some preferred embodiments of the present 
invention, there is provided a method for spatially manipulating polarization- 
20 dependent geometrical phases of an incident light beam, the method comprising: 
providing a substrate comprising a plurality of zones of gratings with discretely 
varying orientation, the orientation denoted by tf(x,y)that is equal to half of a desired 
geometrical phase (DGP) modulus 2tt, each grating with a local period that is 
substantially smaller than the wavelength of the incident light beam; 
25 irradiating the light beam on the substrate. 

Other aspects of the invention will be appreciated and become clear after 
reading the present specification and reviewing the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
30 In order to better understand the present invention, and appreciate its practical 

applications, the following Figures are provided and referenced hereafter. It should be 
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noted that the Figures are given as examples only and in no way limit the scope of the 
invention. Like components are denoted by like reference numerals. 
Figure 1: Illustration of the principle of PBOEs by use of the Poincare sphere. 
Figure 2: The geometry of the space-^variant subwavelength grating , as well as the 
5 Diflfractive Geometrical Phase, DGP, for incident \ R) and | L) polarizations. 

Figure 3: Measurements of the transmitted far-field for the subwavelength grating 
PBOE, when the retardation is ^;zx2and ^irrespectively, for incident (a) circular 
left, (b) circular right and (c) linear polarizations. 

Figure 4: Illustration of polarization diffraction gratings fabricated using 
10 subwavelength quasi-periodic structures. The orientation of the subwavelength 
grooves, 0{x) , varies periodically in the x-direction, resulting in an element with 
effective birefringence, whose optical axes (marked by the dark and light arrows in 
the picture), rotate periodically. The polarization diffraction grating has period d, 
which is larger than the incident wavelength, whereas the local subwavelength 
1 5 period is A < A . The local optical axes at each point are oriented parallel and 
perpendicular to the subwavelength grooves. 

Figure 5: A diagram describing the operation of polarization diffraction gratings. A 
beam with polarization |E,„) is incident on the polarization grating. The resulting 
beam comprises three polarization orders, EPO( E Polarizing Order), which maintains 

20 the original polarization and does not undergo phase modification. The RPO (Right 
Polarizing Order) that is |R) polarized, and whose phase is modified by 20{x) , and 
LPO (Left Polarizing Order) that is |L) polarized, and whose phase is modified by 
- 20ix) . Since ^(jc) is periodic, the RPO and the LPO undergo diffraction, resulting 
in the appearance of discrete diffraction orders. 

25 Figure 6: (a) The magnified geometric representation of the continuous blazed 

polarization diffraction grating; (b) the resulting DGPs for RPO and LPO formed by 
this structure, (c) shows a magnified image of a region on the subwavelength 
structure, demonstrating the local subwavelength period A(x,;;), the local 
subwavelength groove orientation ffix^y) , and the local subwavelength grating vector 

30 Kg . (d) Scanning electron microscope image of the GaAs dielectric structure. 
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Figure 7: Images and experimental (dots) and calculated cross-section (solid curves) 
of the far-field of the beam transmitted through the blazed polarization diffiraction 
gratings, when the incident beam has |L) polarization, when the incident beam has 

|R) polarization, and when the incident beam has linear polarization ^l"!^) * 

S Figure 8: Measurements and predicted diffraction efRciencies in the 1^\ -1^' and 0 
orders of the metal-stripe and dielectric blazed polarization gratings for various 
incident polarizations. The different incident polarization were achieved by rotating a 
quarter waveplate (QWP) in front of the linearly polarized light emitted from the CO2 
laser. The graphs show the efficiencies as a function of the orientation of the QWP. 
10 The efficiencies are normalized relative to the total transmitted intensity for each 
grating. 

Figure 9: geometric representation of the subwavelength structure of the binary 
polarization diffraction grtaing (middle), as well as the DGP that results from this 
structure (top). The pictures (bottom) are Scanning Electron Microscope images of the 
15 dielectric subwavelength structure. 

Figure 10: Measured intensity and polarization in the various diffi-action orders of the 
binary polarization diffraction grating for incident (a) |L) polarization, (b) incident 

|R) polarization and (c) incident j polarization, as well as the intensity transmitted 
through a combination of the binary diffraction grating and a circular polarizer 
20 oriented to transmit |R) polarization for (d) incident |L) polarization (e) incident 

|R) polarization and (f) incident polarization. The intensities are normalized so 
that the maximum intensity in each graph is equal to 1 . 

Figure 11: The far-field images and calculated and measured cross-sections of the 
beam transmitted through the circular symmetric polarization mode switching, (a) 
25 when the incident polarization is |L) , as well as (b) the image and cross-sections of 

the transmitted |L) component, and (c) the image and cross sections of the 

transmitted |R) component. Also shown, (bottom), the geometry of the 

subwavelength quasi-periodic structure as well as the spiral DGP caused by this 
element. 
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Figure 12: (a) Illustration of the principle of PBOEs by use of the Poincare sphere; 
the insets display the different local orientations of subwavelength grooves; (b) 
illustrates the geometrical definition of the grating vector. 
Figure 13: The geometry of the subwavelength gratings for different topological 
5 charges, as well as an image of typical grating profile taken with a scanning-electron 
microscope. 

Figure 14: (a) The interferogram measurements of the spiral PBOEs; (b) the 

corresponding spiral phases for dififerent topological charges. 

Figure 15: Experimental far-field images as well as their calculated and measured 
10 cross sections for the helical beams with / = 1 - 4 

Figure 16: Actual quantized phase Fiifi^) as a function of the desired phase , as 

well as the multilevel discrete binary local grating orientation. Inset, a scanning 

electron microscope image of a region on the subwavelength structure. 

Figure 17: The magnified geometry of the grating for TV = 4 , as well as the predicted 
1 5 geometrical quantized phase distribution. Also shown, the measured (triangles) and 

predicted (dashed curve) diffraction efSciency as a function of the number of discrete 

levels, N . 

Figure 18: Illustration of the magnified geometry of a multilevel-PBOE focussing 
lens with iV = 4, as well as the predicted multilevel geometrical phase. Inset, the 
20 image of the focused spot size as well as the measured (dots), and theoretically 
calculated (solid curve) cross section. 

Figure 19: Schematic presentation of the near-field Fourier transform polarimetiy 
based on a quantized space-variant subwavelength dielectric grating followed by a 
subwavelength metal stripes polarizer. 

25 Figure 20: (a) The magnified geometry of the quantized space-variant dielectric 

grating with a number of quantized levels ^ = 4 . (b) The quantized rotation angle of 
the subwavelength grating as a function of thex coordinate; the inset illustrates the 
local groove orientation, (c) Scanning electron microscopy image of a region on the 
subwavelength structure (d) Scanning electron microscopy image of a cross section of 

30 the subwavelength grooves. 
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Figure 21: The diffraction orders emerging from the QSG; zero order (dashed), first 
order for |R) and |L) polarized beams (solid), and higher order for |R) and 

I L) polarized beams (dash-dot). 

Figure 22: Measured (dotes) and predicted (solid curves) of the normalized 
5 transmitted intensity as a function of the x coordinate along the QSG, when the fast 
axis of the rotating QWP was at angles of (a) 0® , (b) 20° , and (c) 45° ; the insets 
show the experimental images of the of the near-field transmitted intensities. 
Figure 23: Measured (dotes) and predicted (solid curves) values of the normalized 
Stokes parameters, (a) S^/S^, (b) S^/Sq, (c) S^/S^.bssl function of the orientation 
10 of the QWP. 

Figure 24: Measured (dotes) and predicted results (solid curves) for (a) azimuthal 
angle ^ , and (b) ellipticity angle , as a function of the orientation of the QWP. 
Figure 25: The calculated (solid curve) and measured (circles) DOP as a function of 
the intensity ratio of the two independent lasers having orthogonal linear polarization 
15 states, as used in a setup depicted in the top inset. The bottom inset shows calculated 
(solid curves) and measured (dots) intensity cross-sections for the two extremes, Ii=l2 
(DOP=0.059) and l2=0 (DOP=0,975). 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
20 Several embodiments of geometrical phase optical elements with space-variant 

subwavelength gratings and their implementations thereof are hereby presented, 
without derogating generality and without limiting the scope of the present invention 
to these embodiments only. 

25 1. Space-variant Pancharatnam-Berry phase optical elements with sub- 
wavelength gratings: 

Space-variant Pancharatnam-Berry phase optical elements based on computer 
generated subwavelength gratings are presented herein. By continuously controlling 
the local orientation and period of the grating any desired phase element can be 

30 achieved. We present a theoretical analysis, and experimentally demonstrate a 
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Pancharatnam-Beny phase based diffiactioif grating for laser radiation at a 
wavelength of 10.6 //m. 

We consider optical phase elements based on the space-domain Pancharatnam- 
Beny phase. Unlike diffractive and refractive elements, the phase is not introduced 
5 through optical path differences, but results from the geometrical phase that 

accompanies space-variant polarization manipulation. The elements are polarization 
dependent, thereby enabling multi-purpose optical elements, suitable for applications 
such as optical switching, optical interconnects and beam splitting. We show that such 
elements can be realized using continuous computer-generated space-variant 

10 subwavelength dielectric gratings. The continuity of the gratings ensures the 
continuity of the resulting field thereby eliminating diffraction associated with 
discontinuities, and enabling the fabrication of elements with high diffraction 
eflBciency. We experimentally demonstrate Pancharatnam-Berry phase diffraction 
gratings for CO2 laser radiation at a wavelength of 1 0.6fxm^ showing an ability to form 

1 5 complex polarization-dependent continuous phase elements. 

Figure 1 illustrates the concept of Pancharatnam-Berry phase optical elements 
(PBOEs) on the Poincare sphere. Circularly polarized light is incident on a waveplate 
with constant retardation and a continuously space varying fast axis whose orientation 
is denoted by 6(x,y), We will show that since the waveplate is space varying, the beam 

20 at different points traverses different paths on the Poincare sphere, resulting in a 
space-variant phase-front modification originating from the Pancharatnam-Berry 
phase. Our goal is to utilize this space-variant geometrical phase to form novel optical 
elements. 

It is convenient to describe Pancharatnam-Berry Optical Elements (PBOEs) 
25 using Jones calculus. In this formalism, a waveplate with a space-varying fast axis, is 
described by the operator, 

T(jc,;;) = R(^(jc,7))J(^)R"* {9{x,y)) , where is the operator for a waveplate with 
retardation ^, R is the operator for an optical rotator and ^is the local orientation of 
the axis at each point {x,y). For simplicity, we work with the helicity basis in which 
30 1 1) denotes left-hand circular polarization, and \r) denotes right-hand circular 
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polarization. In this representation, Tix,y) has the explicit form, 

T(,.^) = cos(<^/2)[^^ lJ-'^«'(*^/Hexp(-/2^(x.^)) 0 J' 

Thus, if a beam with polarization E,- is incident on Tix^y), the resulting beam has the 

explicit form, 

5 E, =T(x,j.)E, =cos(«i/2)E,-isin(^^/2)[(E,|/J)expH2tf)|i^ 
(2) 

Using Eq. (2) we can calculate the Pancharatnam phase-front of the resulting wave. 
We define this phase-front based on Pancharatnam's definition for the phase between 
two beams of different polarization as 9?^(x,>;) = arg[(Eo(0,0)|Eo(jc,;^))]. This 

10 calculation yields that for incident | R) polarization, 

(x, ^) = + arctan[cos ^> tan ^] = + arctan[sin(2;ir) tan O] , where z is the 
ellipticity of the resulting beam. Geometrical calculations show that tpp is equal to half 
the area of the geodesic triangle, Q on the Poincare sphere defined by the pole | R) , 
|Eo(0)) aiid |Eo(^)) , as illustrated in Fig. 1, yielding the expected Pancharatnam- 

15 Berry phase. Similar results can be found for any incident polarization. 

Consequently, if a circularly polarized beam is incident on a space-variant 
polarization-state manipulator, it is subject to geometrical phase modification. Based 
on Eq. (2), the resulting wave consists of two components, the zero order, and the 
diffracted order. The zero order has the same polarization as the original wavefront, 

20 and does not undergo any phase modification. On the other hand, the diffracted order 
has a polarization orthogonal to that of the incoming wave, and it's phase at each 
point is equal to twice the local orientation of the waveplate, 0(x,y). Since the phase 
modification of the wavefront is purely geometrical in origin, the phase of the 
diffracted orders must also be geometrical. We therefore define the Diffractive 

25 Geometrical Phase (DGP) as the phase of the diffracted orders when the incident 

beam is circularly polarized. For incident \R) and \L) polarizations the DGP is equal 

to •20(x,y) and 26(x,y) respectively. By correctly determining the local orientation of 
the waveplate, any desired DGP can be realized, enabling the realization of phase 
operators such as lenses or diffraction gratings. Furthermore, since the orientation of 
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the waveplate only varies from 0 to ;^ the DGP is defined modulus 27^ and the 
elements are analogous to diffractive optical elements. However, unlike diffractive 
optica] elements the phase modification in PBOEs does not result from optical path 
differences, but from polarization-state manipulation, and is intrinsically polarization 
5 dependent, with the transmission function given by the matrix T(x,y) as defined in Eq. 
(2). This equation also shows that the diffraction efficiency of the PBOEs depends on 
the retardation of the waveplate. 

A case of special interest is ^=n. In this case, we find that the diffraction 
efficiency is 100%, and that | R) polarization is completely converted into 1 1) 

10 polarization. However, despite the fact that the resulting polarization is space- 
invariant. The Pancharatnam phase, = -29{Xyy) , is equal to the DGP. This phase 

corresponds to half the area encompassed by two geodesic paths between the poles 
that form an angle of with one another, as illustrated in Fig. 1 . Thus, the DGP is 
equal to the geometrical Pancharatnam-Berry phase of a PBOE with 100% diffraction 
15 efficiency. Note that PBOEs operate in different ways on the two helic polarizations. 
Consequently, a PBOE lens designed for a wavelength A, with focal length / 
designed by choosing the direction of the waveplate so that 20{x,y) = * 

will be a converging lens for | R) polarization, and a diverging lens for | L) 
polarization. 

20 PBOEs can be realized using space-variant subwavelength gratings. When the 

period of the grating is much smaller than the incident wavelength, the grating acts as 
a uniaxial crystal (Z. Bomzon, G. Biener, V. Kleiner and E. Hasman, ,Opt. Lett. 27, 
1 88 (2002)). Therefore by correctly controlling the depth, structure, and orientation of 
the grating, the desired PBOE can be made. To design a PBOE, we need to ensure that 

25 the direction of the grating stripes, 6{x^y) , is equal to half the desired DGP which we 
denote as qf^ix^y) . Next we define a grating vector 

=i^o(^»3^)[c<^s(^3^(x,>')/2)x + sin(«?rf(jc,j/)/2)y], where xand y are unit vectors 
in the x and y direction, = 2;t/ A(jc,>;) is the spatial frequency of the grating {A is 
the local subwavelength period) and (x, y)/2 is the space- variant direction of the 

30 vector defined so that it is perpendicular to the grating stripes at each point. Next, to 
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ensure the continuity of the grating thereby ensuring the continuity of the resulting 
optical field, we require V x = 0 , resulting in a differential equation that can be 

solved to yield the local grating period. The grating function ^ ^defined so that 
V^^ = ) is then found by integrating Kg over an arbitrary path (Z. Bomzon, V. 

5 Kleiner and E. Hasman, Opt. Lett. 26, 33 (2001)). 

We designed a PBOE that acts as a diffraction grating by requiring that 
(p^ - {jLnld)j\^^^ , where d is the period of the structure. Applying this to the 
grating vector and solving the equation V x = 0 yields, 

=(2;r/Ao)exp-(;rv/rf)[cos(;Dc/rf)x-sin(;cc/rf)y], where Aq is the 

10 subwavelength period aty=0. The grating function is then found as 

^^{x^y) = (2t// Ao)sin(;cc/rf)exp{-;zy/<i) . We then realized a Lee-type binary 

grating describing the grating function, Z. Bomzon, V. Kleiner and E. Hasman, 
Opt. Commun. 192, 169 (2001)). The grating was fabricated for CO2 laser radiation 
with a wavelength of 10.6[im, with A^ = 2//w and rf=2.5/nm, and consisted of 12 

1 5 periods of d. We formed the grating with a maximum local subwavelength period of 
/l=3.2//m because the Wood anomaly occurs at 3.24//m for GoAs. The grating was 
realized onto a 500 thick GoAs wafer using contact printing photolithography and 
electron-cyclotron resonance etching with BCh to nominal depth of 2,5 //m, to yield a 
retardation of ^ = ;r/2 . By combining two such gratings we obtained a grating with 

20 retardation ^ - ;r . Figure 2 illustrates the geometry of the grating, as well the DGP 

for incident |£) and \B!) polarization states as calculated firom Eq. (2). The DGPs 

resemble blazed gratings with opposite blazed directions for incident | and | 

polarization states, as expected from our previous discussions. 

Following the fabrication, we illuminated the PBOEs with circular and linear 
25 polarization. Figure 3 shows the experimental images of the diffracted fields for the 
resulting beams, as well as their cross-sections when the retardation was ^=^and 
^==;r respectively. When the incident polarization is circular, and ip=7c/2^ close to 
50% of the light is diffracted according to the DGP into the first order (the direction 
of diffraction depends on the incident polarization), whilst the other 50% remains 
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undiffracted in the zero order as expected from Eq. (2). The polarization of the 
diffracted order has switched helicity as expected. For ^=7i^ no energy appears in the 
zero order, and the diffraction efficiency is close to 100%. When the incident 
polarization is linear E,. = ^ (j ^) + 1 1)), the two helic components of the beam are 

5 subject to different DGPs of opposite sign, and are diffracted to the first orders in 

different directions. When ^^=^2^ the zero order maintains the original polarization, in 
agreement with Eq. (2), whereas for retardation n the diffraction is 100% efficient for 
both circular polarizations, and no energy is observable in the zero order. 

To conclude, we have demonstrated novel polarization-dependent optical 

10 elements based on the Pancharatnam-Berry phase. Unlike conventional elements, 
PBOEs are not based on optical path difference, but on geometrical phase 
modification resulting from space-variant polarization manipulation. 



15 

2. Polarization beam-splitters and optical switches based on space-variant 
computer generated subwavelength quasi-periodic structures: 

Polarization beam-splitters and optical switches based on subwavelength 
quasi-periodic structures are presented. By locally controlling the orientation and 

20 period of the subwavelength grooves, birefringent elements for which the optical axes 
vary periodically, are realized. We present a theoretical discussion of these elements, 
as well as a detailed description of the design and realization procedures. We show 
experimental results for infra-red radiation at a wavelength of 10.6|im. 

Polarizing beam-splitters are essential components in polarization-based 

2S systems such as ellipsometers, magneto-optic data storage and polarization-based 
light modulators. Often these applications require that the elements provide high 
extinction ratios over a wide angular bandwidth while maintaining compact and 
efficient packaging. Conventional polarizing beam-splitters, employing either natural 
crystal birefringence or polarization-sensitive multilayer structures are usually, either 

30 cumbersome or sensitive to angular change, and therefore do not fully meet these 
requirements. 
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Contemporary research has begun to address the use of polarization diffraction 
gratings as beam-splitters and optical switches. Unlike scalar diffraction gratings that 
are based on periodic modification of phase and amplitude, polarization diffraction 
gratings introduce a periodic spatial change of the state of polarization leading to 
5 polarization-dependent diffraction. Furthermore, the polarization of the diffracted 
orders is generally different from that of the incident beam. Such a device was 
demonstrated by Davis et al, who used liquid crystals to create a waveplate with 
space-varying retardation (J. A. Davis, J. Adachi, C. R. Femtodez-Pousa and 1. 
Moreno, Opt. Lett. 26, (2001) 587-589). Alternatively, Gori suggested a grating 

10 consisting of a polarizer with a spatially rotating transmission axis (F. Gori, Opt. Lett. 
24, (1999) 584-586) and Tervo and Turunen suggested that beam-splitters consisting 
of spatially rotating wave-plates could be realized using subwavelength gratings (J. 
Tervo and J. Turunen, Opt Lett. 25, (2000) 785-786). More recently a polarization 
diffraction grating based on spatially rotating nematic liquid crystals has been 

15 demonstrated (B. Wen, R. G. Petschek and C. Rosenblatt, Appl. Opt. 41, (2002) 
1246-1250). 

Space- variant polarization-state manipulations using computer-generated 
subwavelength structures were disclosed (Z. Bomzon, V. Kleiner and E. Hasman, 
Opt. Commun 192, (2001) 169-181, Z. Bomzon, V. Kleiner and E. Hasman, Opt. 

20 Lett. 26, (2001) 33-35, and Z. Bomzon, G. Biener, V. Kleiner, E. Hasman, Optics 
Lett. 27, (2002) 285-287). When the period of a subwavelength periodic structure is 
smaller than the incident wavelength, only the zeroth order is a propagating order, and 
all other orders are evanescent. The subwavelength periodic structure behaves as a 
uniaxial crystal with the optical axes parallel and perpendicular to the subwavelength 

25 grooves. Therefore, by fabricating quasi-periodic subwavelength structures, for which 
the period and orientation of the subwavelength grooves was space- varying, we 
realized continuously rotating waveplates and polarizers for CO2 laser radiation at a 
wavelength of 10.6^m. Furthermore, we showed that such polarization manipulations 
necessarily led to phase modification of geometrical origin, which left a clear 

30 signature on the propagation of the resulting wave. The phase introduced did not 

result from optical path differences but solely from local changes in polarization and 
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was in-&ct a manifestation of the geometrical Pancharatnam-Beny phase (Z. 
Bomzon, V. Kleiner and E. Hasman, Opt. Lett. 26, (2001) 1424-1426). 

Disclosed herein is an experimental demonstration of polarization diffraction 
gratings based on space-variant computer-generated subwavelength structures. We 
5 present a novel interpretation of these elements and show that the polarization-related 
diffraction is indeed connected to the space-vaiying Pancharatnam-Berry phase 
mentioned above. We present experimental results for infira-red CO2 laser radiation 
that include a polarization diffraction grating based on a space-variant continuous 
metal-stripe subwavelength structure, a continuously rotating dielectric 

10 subwavelength structure and a binary waveplate, for which the direction of the 

subwavelength grooves varies discretely. We also demonstrate a circular symmetric 
polarization mode switching based on a computer-generated subwavelength structure, 
thus enabling alternation between intensity distributions of a bright and daric center. 
Figure 4 is a schematic representation of the use of a subwavelength structure 

IS as a polarization diffraction grating. The local orientation of the subwavelength 
grooves, 0{x) , varies linearly in the x-direction, to form a polarization diffraction 
grating comprising a birefringent element with optical axes (which are parallel and 
perpendicular to the grating grooves) that rotate periodically in the x-direction. The 
period of the polarization diffraction grating, is larger than the incident wavelength, 

20 A ,whereas the local subwavelength period of the grooves, Aix^y) , is smaller than the 
incident wavelength. When a plane-wave with uniform polarization is incident on 
such a periodic subwavelength stmcture, the transmitted field will be periodic in both 
polarization and phase, therefore, we can expect this field to yield discrete diffraction 
orders in the far-field. 

25 It is convenient to describe subwavelength quasi-periodic structures such as 

the one depicted in Fig. 4 using Jones Calculus. In this representation, a uniform 
periodic subwavelength structure the grooves of which are oriented along the y-axis 
can be described by the Jones matrix, 

^=(0 • <■> 

30 where , are the real amplitude transmission coefficients for light polarized 

perpendicular and parallel to the optical axes and ^ is the retardation of the grating. 
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If the orientation of the subwavelength grooves is space-varying, i.e. different at each 
location, then the subwavelength structure can be described by the space-dependent 
matrix, 

Tc(x)^M{a(x))jM^\0(x)) , (2) 

f COS0 -sin^^ . 

where 9(x) is the local orientation of the optical axis and M(^) = ^ \^^^ 

\^sin0 COS0 J 

two-dimensional rotation matrix . 

For convenience we adopt the Dirac bra-ket notation, and convert T^(x) to 



the helicity base in which | R) = j and | L) = j are the two-dimensional 



unit 



vectors for right-hand and left- hand circularly polarized light. In this base, the space- 
10 variant polarization operator is described by the matrix, T(x) = 11X^11"' ,where 

U = ^ M is a unitary conversion matrix. Explicit calculation ofT(x) yields, 

Thus for an incident plane- wave with arbitrary polarization |E„, ) we find that the 
resulting field is, 
15 |E^,> = V^|E,„>H.V;^e«-<->|R> + V^e--^->]^ , (4) 

where tIe =^\W.+ty^'*i[ > -/,e'*)(E,„ |L)f ,7, =|i(^ - v'')(E,>, M'. 

are the polarization order coupling efficiencies and (^|>9) denotes inner product. 

Figure 5 is a graphic representation of the results of Eq. (4). It shows that 
I E^„, ) comprises three polarization orders ; the | E;„ ) polarization order (EPO), the 

20 I R) polarization order (RPO) and the | L) polarization order (LPO). The EPO 

maintains the polarization and phase of the incident beam, whereas the phase of the 
RPO is equal to 2d{x) , and the phase of the LPO is equal to - 20{x) . We note that 

the phase modification of the |R)and the |L) polarization orders results solely from 

local changes in polarization and is therefore geometrical in nature. We therefore 
25 denote this phase as the Diffractive Geometrical Phase (DGP) (Z. Bomzon, G. Btener, 
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V. Kleiner and E. Hasman, "Space-variant Pancharatnam Berry phase optical 
elements with computer-generated subwavelength gratings", Opt. Lett, in press, 
(2002)). 

The DGP for the | R) polarization order is opposite in sign to that of the |L) 

5 polarization order, and if ^(jc) is periodic, the functions e'^^^'^ and e''^^^'^' that appear 
in Eq. (4) can be developed into Fourier series. Taking into account the connection 
between the Fourier series of a function and the Fourier series of its complex 
conjugate, this leads to the equation, 

\E^.)-^K)^4^t^^^^ (5) 

10 where a„ = Je''^^'>e''''""'''dx . 

Thus we find that the diffraction efficiency into the m ^ order of the RPO 

(17^ = \a„ 1^ ), is equal to the diffraction efficiency into the ^n * order of the LPO 

{Tf!:^ = |«« f )y conclude that the RPO and LPO are diffracted in opposite senses. 

Based on Eqs. (4-5) we find that there are three degrees of freedom associated 
15 with the design of polarization diffraction gratings. The first degree of fi-eedom is the 
determination of the subwavelength structure parameters, and ^ . These 

parameters determine the amount of energy coupled into the EPO. The second degree 
of freedom is the grating orientation 0{x) . 0(x) that determines the DGP, thereby 
determining the diffraction efficiency into all diffraction orders. The third degree of 
20 freedom is the incident polarization |E^„ ) . It determines the ratio between the energy 

in the RPO and the energy in the LPO. In the next paragraphs, we intend to 
demonstrate how these three degrees of freedom can be utilized for the design and 
realization of polarization beam-splitters and optical switches using subwavelength 
quasi-periodic structures. 
25 Supposing we wish to design a blazed polarization diffraction grating, i.e. a 

grating for which all the diffracted energy is in the 1*' order when the incident beam is 
|R) polarized. For |E,„) = |R) we find that fJx=0. Consequently the transmitted 

beam only consists of the LPO and the EPO. Since the EPO does not undergo any 
phase modification, all of its energy is located in the zero order, and the only order 
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that contributes energy to the order is the LPO. In order to ensure that all the 
energy of the LPO will be diffracted into the I'** order, it is required that the DGP for 
the LPO be equal to 2m /d\ • Consequently, we find that 0ix) = -m/d\ . 
Next, in order to ensure that no energy is found in the zero-order, we require that 
5 7£ = 0 . This condition leads to the solution = ty and ^ = 7r . Thus by determining 

the incident polarization, |E;„) , the grating orientation, 0{x) , and the grating 
parameters t^^ = and ^ = ;r , we are able to create the desired diffraction pattern. 
In addition, we note that for such a grating, the DGP for the RPO is 
20(x) = -2ml d\ , and therefore if |E^) = |L), the grating is blazed in the 

10 opposite direction. Thus for arbitrary incident polarization, the diffracted energy will 
be distributed between the 1** and -1*' orders. The distribution is dependent on the 
polarization of the incident beam. Furthermore, the polarization in the 1*^ order will 
always be | L) , and the polarization in the -1 ^ order will always be | R) . Thus, by 
switching the incident polarization between an | L) state and an | R) state, an optical 

1 5 switch can be realized. Furthermore, if we choose (p^n ^ then some of the incident 
energy will be coupled into the EPO, resulting in the appearance of a zero-order that 
maintains the polarization of the incident beam, thereby demonstrating the usefulness 
of such a device as a variable polarization-dependent beam splitter. 

We now focus our attention on the design of the blazed grating discussed 

20 above using a quasi-periodic subwavelength structure. Since the determination of the 
grating parameters t^, and ^ depend mainly on the subwavelength groove profile, 

and not on the subwavelength groove orientation, we begin by determining the desired 
subwavelength groove orientation, 0{x) , and period, A(x,y) , for the subwavelength 
structure with the desired DGP. The grating parameters t^^t^ and ^ are later 

25 determined by choosing a fabrication process that yields a grating profile with the 
desired birefringence. 

To design a continuous subwavelength structure with the desired DGP, we 
define a subwavelength grating vector, {x^y) , oriented perpendicular to the 

desired subwavelength grooves. 
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K^ix,y) = K^{x,y)cos(fDc/d)x-K^(x,y)s\n(;ix/d)y . (6) 
^oi^>y) = 2;r/ A(jc,>;) is yet to be determined as the local spatial frequency of the 
subwavelength structure. Figure 6(c) illustrates this definition of {x^y) • To ensure 
the continuity of the subwavelength grooves, we require that V x = 0 , leading to 
S the partial differential equation, 

— ^cos(;Dc/i/) + — 2-sin(;cc/£/) + ^A:oCos(;Dc/f/) = 0 , (7) 
cy dx * ^ 

with the boundary condition K^ixSS) = 2;r/ , where Aq is the local subwavelength 

period at>^. The solution to this problem is given by, 

K^(jc,y) = -^exp(-;zy/rf)[cos(;zx/rf)x -sin(;cc/rf)y] . (8) 

1 0 Consequently, the grating function is then found by integrating Kg over an arbitrary 
path to yield, 

'id 

' ^gix,y) = —sinim/d)exp('-;ry/d)(9) 
Ao 

We realized a Lee-type binary subwavelength structure mask (see W. H. Lee, 
Appl. Optics 13, (1974) 1677-1682) described by the grating function of Eq. (9) using 
1 5 high-resolution laser lithography. The amplitude transmission for such a Lee-type 
binary mask can be derived as, 

t{x,y) = U,[cosiM - cos(;r ^)] , (10) 
where is the unit step function defined by, 

' ' 1 0, 7<0 

20 and where q is the duty cycle of the grating which was chosen as 0.5. Figure 6 

illustrates the geometry of a Lee-type binary subwavelength, as well as the resulting 
DGPs for the RPO and the LPO formed by this structure. The figure shows a 
continuous quasi-periodic subwavelength structure with a local subwavelength period 
A(x, J/) where at each location on the element, the grooves are oriented perpendicular 

25 to the required fast-axis, resulting in the desired polarization diffiaction grating. The 
resulting DGPs resemble the phase function of a scalar blazed grating. The RPO is 
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blazed in the opposite direction of the LPO as discussed above. Hence, incident | R) 
polarization is diffracted in the opposite direction of incident |L) polarization. The 

local subwavelength periodicity gives the structure its birefringence, whilst the 
continuity of the subwavelength grooves ensures the continuity of the resulting field. 
5 Furthermore, we note the space-varying nature of A(jc, y) , This is a necessary result 
for the requirement of continuity posed on flie subwavelength grooves. 

We realized three different subwavelength structures with the geometry shown 
in Figure 6. The first element was realized as a metal stripe subwavelength structure 
using contact photolithography and lift-off, and the other two elements were dielectric 

10 gratings realized using contact photolithography and dry etching techniques. The 

elements were realized for CO2 laser radiation at a wavelength of \0.6\im on 500 nm 
thick GaAs wafers. We fabricated the gratings with = , and d = 2,5mm , 
consisting of 12 periods of d . The gratings were formed with maximum local 
subwavelength period of A = 3.2 /im because the Wood anomaly occurs at 

IS 3.24 jm for GaAs. The metal stripes consisted of a lOnm adhesion layer of Ti and 

60nm Au with a duty cycle of 0.6 yielding measured values of = 0.6 , = 0.2 , and 
^ = 0.6;r . The dielectric gratings were fabricated using electron-cyclotron resonance 
etching with BCI3 to nominal depth of 2.5 //wi and duty cycle of 0.5, resulting in 
measured values of a retardation ^ = ;t/2 , and = = 0.9 . By combining two such 

20 gratings, we obtained a grating with retardation <f>~7t, and =t^ = 0.89 . These 

values are close to the theoretical predictions achieved using rigorous coupled wave 
analysis (M.G. Moharam and T.K. Gaylord, J. Opt. Soc. Am. A 3, (1986) 1780-1787). 
Figure 6(d) shows a scanning electron microscope image of one of the dielectric 
structures. We note the local periodicity of the structure and the clear profile of the 
25 subwavelength grooves. 

Following the fabrication we illuminated the structures with polarized light. 
Figure 7 shows images of the far-field intensity distributions of the transmitted 
beams, as well the measured and predicted intensity cross-sections for incident |R) 

polarization, for incident |L) polarization and for incident linear polarization 
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( I ) = 1 ). There is a good agreement between experiment and theory, which was 
calculated using Eq. (4), and far-field Fraunhoffer intensities. For incident |R) 

polarization only the zero order and the 1** order appear. The zero order is due to the 
EPO, and maintains the polarization of the incident beam, whereas the 1^' order is 
5 derived from the LPO and has |L) polarization. For incident |L) polarization, the 

zero-order and the -1*' order appear, and the polarization of the -1^ order is | R) . In 
the case of incident linear polarization, we note the linear polarization of the zero 
order, the | L) polarization of the 1 order and the | R) polarization of the -P^ order as 
discussed above. We note that for the grating in which t^^ = t^, = 0.89 and (^ = 7ty 
10 7^ = 0 , consequently, all of the energy is diffracted into the 1** and -1 orders, as 
expected. 

Figure 8 shows the predicted and measured diffraction efficiencies for the 
three blazed polarization diffraction gratings for various incident polarization states. 
The diffraction efficiencies are normalized relative to the total transmitted intensity. 

1 5 The different polarization-states were achieved by rotating a quarter wave-plate in 

front of the linearly polarized light emitted from the laser. The experiments agree with 
the predictions. The diffraction efficiency in the zero-order is equal to Tfg, It has a 
different value for each of the grating, however, for each grating, it remains constant 
regardless of the incident polarization. On the other hand, the diffraction into the 1** 

20 and -1^' orders depends on the incident polarization, illustrating the usefulness of 
polarization diffraction gratings as variable polarization beam-splitters and light 
modulators. 

To further examine the use of polarization diffraction gratings as beam 
splitters and optical switches, we fabricated a binary polarization diffraction grating 
25 using a subwavelength dielectric structure. Each period, (d >XX of the grating is 
comprised of two regions. The subwavelength grating vector in the first region 
pointed along the x -axis, and the subwavelength grating vector in the second region 
pointed along the y -axis, i.e, 

^ ' \til2 dl2<x<d ^ ' 
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We fabricated the subwavelength structure on a GaAs wafer with = 0.95 , 
ty = 0.84 , and ^ = 0.45;r . The structure was fabricated with d=200[im and 
A(x^y) = 2//m . Figure 9 shows a schematic representation of the subwavelength 
structure, a graph depicting the resulting DGP and scarming electron microscope 
images of the subwavelength structure which we had fabricated. We note that the 
DGP is the same for both the RPO and the LPO. It resembles a scalar binary tt -phase 
grating. Consequently the diffraction efficiencies for both the RPO and the LPO, will 
be the same as those of a scalar binary ;r -phase grating (E. Hasman, N. Davidson, 
and A.A. Friesem, Opt. Lett. 16, (1991) 423-425), i.e.. 



10 ^-=-.- = ""^7^^^ . (12) 



Note that for a binary ;r -phase grating only the odd orders appear and consequently 

=tJq =^0. Furthermore, Eq. (12) yields that fj{"^ = T}t^ = 0.405 , and 
^L,* =7.Y =0.045. 

Figure 10 shows the measured diffraction efKciencies when the incident beam 
15 has (Figs. 7(a,b)) circular and (Fig. 10(c)) linear polarization, as well as the efficiency 
when the transmitted beam has passed through a circular polarizer oriented to transmit 
|R) polarized light (Figs. 10(d-f)). We note that when the beam does not pass 

through a circular polarizer, the intensity in the various diffracted orders is not 
dependent upon the incident polarization, however the polarization-state of the 

20 diffracted orders does depend on the incident polarization. We note that the intensity 
of the diffracted orders on the right is equal to the intensity of the diffracted orders on 
the left, indicative of synunetrical phase structure of the DGP. Furthermore the ratio 
between the intensity in the 1** and 3"* orders is 7| /773 « 9 , in agreement with Eq. 
(12), providing further verification of the binary tt -phase of the DGP. In addition, for 

25 incident | R) polarization, the experimental ratio between the intensities in the zero 
order and in the 1^* order is 3.376. This agrees with the predicted ratio of 

^£ linuVx ) = [ I ^ + f / 0-405| f , - f ] . as predicted using Eqs. (4, 1 2). 
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When a circular polarizer, oriented to transmit only | R) polarized light is applied to 
the beam, we notice that for incident |R) polarization (Fig* 10(e)) only the zero order 
appears, (this is because rj^ -0). For incident |L) polarization (Fig. 10(d)) only the 
orders other than the zero order appear (this is because the EPO has |L) polarization), 

5 and for incident linear polarization (Fig. 10(f)) all diffracted orders appear (this is due 
all three orders being linearly polarized). Thus, by placing a polarization modulator 
such as a liquid crystal cell in front of a setup, comprising a polarization-diffraction 
grating and a circular polarizer, an optical switch could be assembled for applications 
such as optical interconnections in communications. 

10 Until now we have discussed polarization beam splitting and optical switching 

by use of polarization diffraction gratings. However, sometimes a different geometry 
is required. Suppose for instance, we wish to create an optical switch that enables 
switching between an optical circularly symmetric mode with a bright center, and an 
optical circularly symmetric mode with a dark center. This can be done with a quasi- 

1 5 periodic subwavelength structure for which, ^(x,y) = eo(x, j') + c , where 

a> = arctan(;;/x) is the azimuthal angle, and c is a constant number. The DGP of the 
RPO for such an element is equal to 2(a>ix, y) + c) . Thus the RPO carries a vortex 
with a topological charge of 2 (D. Rozas, C. T. Law and G. A. Swartzlander Jr., J. 
Opt. Soc, Am. B. 14, (1997) 3054-3065), and therefore it has a dark center. 

20 Furthermore, since the EPO does not undergo any phase modification, its topological 
charge is zero, and it exhibits a bright center. Therefore, if we design a quasi-periodic 
subwavelength structure with 0ix,y) = fi>(x,>;) + c, choose ^ ;r , and illuminate it 
with |l) polarization, the resulting beam will comprise an |R) polarized vortex 
carrying beam with a dark center and an |L) polarized beam with a bright center. We 

25 can switch between the two modes using a circular polarizer. 

We realized such a quasi-periodic subwavelength structure on a GaAs wafer, 
with 0(x,y) = o)(x^y) + ;r /4, yielding a grating function 

= (2;zro/A/2 Ao)[ln(r/ro) where r is a radial coordinate. We chose 

Aq = 2/4m and = Smm , so that 5mm < r < 8mm and 2jum < A < 3.2/ot, and 
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fabricated a dielectric grating with retardation of ^ = ;r/2 . Note that for incident |l) 
polarization, the transmitted beam had radial polarization in the near field, and that for 
incident |R) polarization, the near field had azimuthal polarization. Figure 11 shows 
far-field images and the measured and calculated cross-sections of (a) the transmitted 
S beam when the incident polarization is |L) as well as the far-field image of the (b) 

|l) and (c) |R) components of the beam as obtained with a circular polarizer. Figure 

11 also shows the geometry of the subwavelength quasi-periodic structure as well as 
the DGP caused by this element. We note the clear vortex in the DGP. The dark 
center of the vortex is evident in the measured results where we note that the |l) 
10 component of the transmitted beam (the EPO), has a bright center without undergoing 
any phase modification, whereas the |r) component (RPO) has a dark center, clearly 

indicative of its topological charge. The results clearly demonstrate the circularly 
symmetric mode switching, which can be realized using subwavelength periodic 
structures. 

15 To conclude, we have demonstrated polarization diffraction gratings as 

polarization-sensitive beam-splitters, as well as optical switches. We have 
demonstrated that the application of subwavelength quasi-periodic structures, for this 
purpose, is not limited to optical switches based on linear polarization diffraction 
gratings, and that more complex designs are possible. The introduction of space- 

20 varying geometrical phases through space-variant polarization manipulations, enables 
new approaches for beam-splitting and the fabrication of novel polarization-sensitive 
optical elements. 

3« Formation of iieiicai beams using Pancharatnam-Berry piiase optics elements: 
25 Spiral phase elements with different topological charges based on space- 

variant Pancharatnam-Beny phase optical elements are presented. Such elements can 
be realized using continuous computer-generated space-variant subwavelength 
dielectric gratings. We present a theoretical analysis and experimentally demonstrate 
spiral geometrical phases for infra-red radiation at a wavelength of 10.6// m . 
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Recent years have witnessed a growing interest in helical beams that are 
exploited in a variety of applications. These include trapping of atoms and 
macroscopic particles, transferring of orbital angular momentum to macroscopic 
objects, rotational frequency shifting, the study of optical vortices, and specialized 
5 alignment schemes. Beams with helical (or spiral) wavefronts are described by 

complex amplitudes M(r,€»,) oc exp(-//fi>) , where r,€y, are the cylindrical coordinates, 
the radial coordinate and azimuthal angle respectively, and / is the topological charge 
of the beams. At the center, the phase has a screw dislocation, also called a phase 
singularity, or optical vortex. Typically, helical beams are formed by manipulating the 

10 light after it emerges from a laser by superposition of two orthogonal (non-helical) 
beams, or by transforming Gaussian beams into helical beams by means of computer- 
generated holograms, cylindrical lenses and spiral phase elements (SPEs). 

Altematively, a helical beam can be generated inside a laser cavity by 
inserting SPEs into the laser cavity. The common approaches of forming SPE are as 

1 5 refractive or difl&active optical elements using a milling tool, single-stage etching 

process with a gray-scale mask, or with multistage etching process. Such helical beam 
formations are generally either cimibersome or sufiTer from complicated realization, 
high aberrations, low efficiency, or large and imstable setup. 

In the present invention novel spiral phase elements (SPEs) are disclosed 

20 based on the space-domain Pancharatnam-Berry phase. Unlike diflfractive and 

refractive elements, the phase is not introduced through optical path differences, but 
results from the geometrical phase that accompanies space-variant polarization 
manipulation. We show that such elements can be realized using continuous 
computer-generated space-variant subwavelength dielectric gratings. Moreover, we 

25 experimentally demonstrate SPEs with different topological charges based on a 
Pancharatnam-Berry phase manipulation, with an axial symmetric local 
subwavelength groove orientation, for CO2 laser radiation at a wavelength of 
10.6 //m. 

The formation of PBOE with spiral geometrical phase, indicates an ability to 
30 form complex continuous PBOEs. 

The PBOEs are considered as waveplates with constant retardation and 
continuously space varying fast axes, the orientation of which is denoted by 0(r^cd). 
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It is convenient to describe PBOEs using Jones calculus. The space-dependent 
transmission matrix for the PBOE is given by applying the optical rotator matrix on 
the Jones matrix of the subwavelength dielectric grating to >deld in helical basis, 

5 where , ty are the real amplitude transmission coefficients for light polarized 

perpendicular and parallel to the optical axes, and is the retardation of the grating. 
Thus for an incident plane-wave with arbitrary polarization |E;„) we find that the 
resulting field is, 

lE„.) = V^|E,„) + ^e-^<-^|R) + ^^^^ , (2) 

10 where rj, = +/,^'^)|', rj, = ~t,.e''){E,„ \hf ,tj, =|K^. -/,^''**)(E„ |R)|', 
are the polarization order coupling efficiencies , denotes inner product, and 

I R) = (l O)'^ and | L) = (O 1^ represent the right-hand and left-hand circular 

polarization components respectively. From Eq. (2) one can see that the emerging 
beam from a PBOE comprises three polarization orders. The first maintains the 
1 5 original polarization state and phase of the incident beam, the second is right-hand 
circular polarized and has a phase modification of 20{r^a)) , and the third is opposite 
to the second at polarization hand and at phase modification. Note that the 
polarization efficiencies depend on Ae groove shape and material as well as on the 
polarization state of the incident beam. For the substantial case of t^-t^-l, and 

20 ^ = ;r an incident wave with | R) polarization is subject to entire polarization state 
conversion and results in emerging field, 

|E<«) = e-«^<-^|L). (3) 

An important feature of Eq. (3) is that the phase factor depends on the local 
orientation of the subwavelength grating. This dependence is geometrical in nature 
25 and originates solely fi-om local changes in the polarization state of the emerging 
beam. This can be illustrated using the Poincare sphere, representing a polarization 
state by three Stokes parameters Si, S2 and S3 as depicted in Fig.l2(a). The incident 
right polarized and the transmitted left polarized waves correspond to the north and 
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the south poles of the sphere, respectively. Since the subwavelength grating is space 
varying, the beam at different points traverses different paths on the Poincare sphere. 
For instance, the geodesic lines A and B represent different paths for two waves 
transmitted through element domains of different local orientations ff(rfi) and 
5 0(ryCo) , Geometrical calculations show that the phase difference 

=-2^(r,fi>)between states corresponding to ^(r,0) and 6{r^a>) orientations, is 

equal to half of the area Q enclosed by the geodesic lines A and B . This fact is in 
compliance with the well-known rule, proposed by Pancharatnam, for comparing the 
phase of two light beams with different polarization, and can be considered as an 
10 extension of the mentioned rule into the space-domain. 

In order to design a continuous subwavelength structure with the desired phase 
modification, we define a space-variant subwavelength grating vector K^(r,a>) , 

oriented perpendicular to the desired subwavelength grooves. Figure 12(b) illustrates 
this geometrical definition of the grating vector. To design a PBOE with spiral 
15 geometrical phase, we need to ensure that the direction of the grating grooves is given 
by 0{r^oji) = /fi?/2 , where / is the topological charge. Therefore, using Fig, 12(b), the 
grating vector is given by, 

K^(r,fi>) = A:o(r,fii){cos[(//2-l)a>]r+sin[(//2-l)a>]fii} , where 

= ItcI h{r^6)) is the local spatial fi-equency of grating with local period of 

20 A(r, fi>) . To ensure the continuity of the subwavelength grooves, we require that 
V X = 0 , resulting in a differential equation that can be solved to yield the local 
grating period. The solution to this problem yields K^ir) - (Itt/ A^Xr^/ rY^^ , where 
Aq is the local subwavelength period at r = . Consequently, the grating function 
(defined so that = ^^g) *s then found by integrating K^(r,£y) over an 

25 arbitrary path to yield, 

^^(r,fi)) = (2;rro/Ao)(ro/r)'''-'cos[(//2-l)fi;]/[//2-l] for 1^2, md 
<t>^(r^co) = (2^r^j / Ao)ln(r/ Tq) for / = 2 . We then realized a Lee-type binary grating 
describing the grating function, , for / = 1,2,3,4 . The grating was fabricated for CO2 
laser radiation with a wavelength of 1 0.6 //m, with = 2//m , /q = 4.7/wm, and a 
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maximum radius of 6mm resulting in 2/im < A(r) < 3.2///n . We formed the grating 
with a maximum local period of 3.2// m in order not to exceed the Wood anomaly of 
GaAs. The magnified geometry of the gratings for different topological charges are 
presented in Fig. 2. The elements were fabricated on 500// m thick GaAs wafers using 
5 contact photolithography and electron-cyclotron resonance etching with BCI3 to 

nominal depth of 2.5//m , resulting in measured values of retardation ^ = ;r /2 , and 
z^ty = 0.9 .These values are close to the theoretical predictions achieved using 

rigorous coupled wave analysis. The inset in Fig. 2 shows a scanning electron 
microscope image of one of the dielectric structures. 
10 Following the fabrication, the spiral PBOEs were illuminated with a right hand 

circularly polarized beam, |R) , at 10.6//m wavelength. In order to provide 

experimental evidence of the resulting spiral phase modification of our PBOEs, we 
used "self-interferogram" measurement using PBOEs with retardation ^ = ;r/2 . For 

such elements the transmitted beam comprises two different polarization orders; | R) 
15 poloarization state, and |l) with phase modification of - according to Eq. (2). 

The near-field intensity distributions of the transmitted beams followed by linear 
polarizer were then measured. Figure 13(a) shows the interferogram patterns for 
various spiral PBOEs, The intensity dependence on the azimuthal angle is of the form 
/ oc 1 + cos(/c£>) , whereas the number of the fiinges is equal to / , the topological 

20 charge of the beam. Figure 13(b) illustrates the phase fi-onts resulting from the 

interferometer analysis, indicating spiral phases with different topological charges. 

Figure 15 shows the far-field images of the transmitted beams having various 
topological charges, as well as the measured and theoretically calculated cross 
sections. We achieved the experimental result by focusing the beam through a lens 

25 with a 500-nun focal length followed by a circular polarizer. We used the circular 
polarizer in order to transmit only the desired |l) state, and to eliminate the 

I R) polarization order that appeared because of the error in the etched depth of the 
grating to achieve the desired retardation of ^ = ;r . Dark spots can be observed at the 
center of the far-field images, providing evidence the phase singularity in the center of 
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the helical beams. We found an excellent agreement between theory and experiment, 
clearly indicating the spiral phases of the beams with different topological charges. 

To conclude, we have demonstrated the formation of helical beams using 
space- variant Pancharatnam-Berry phase optical elements based on computer 
5 generated subwavelength dielectric grating. The formation of the spiral phase by 
PBOE is subject to the control of the local orientation of the grating. This can be 
performed with a high level of accuracy using an advanced photolithographic process. 
In contrast, the formation of SPE based on refractive optics, where the phase is 
influenced by fabrication errors due to inaccuracy of the etched three-dimensional 
1 0 profile. We are currently investigating a photolithographic process to achieve accurate 
control of the retardation phase to yield only the desired polarization order. 

4. Formation of multilevel Pancharatnam-Berry phase diffractive optics using 
space-variant subwavelength gratings: 

1 5 Multilevel discrete Pancharatnam-Berry phase diffractive optics using 

computer-generated space-variant subwavelength dielectric grating is presented. The 
formation of a multilevel geometrical phase is done by discrete orientation of the local 
subwavelength grating. We discuss a theoretical analysis and experimentally 
demonstrate a quantized geometrical blazed phase of polarization diffraction grating, 

20 as well as polarization dependent focusing lens for infra-red radiation at wavelength 
10.6///M. 

One of the most successful and viable outgrowths of holography involves 
diffractive optical elements (DOE). DOEs diffract light from a generalized grating 
structure having nonuniform groove spacing. They can be formed as a thin optical 

25 element that provides unique functions and configurations. High diffraction 

efficiencies for DOEs can be obtained with kinoforms that are constructed as surface 
relief gratings on some substrate. However, in order to reach a high efficiency, it is 
necessary to resort to complex fabrication processes that can provide the required 
accuracies for controlling the graded shape and depth of the surface grooves. 

30 Specifically, in one process a single photomask with variable optical density is 

exploited for controlling the etching rate of the substrate to form the desired graded 
relief gratings, or using multiple binary photomasks so the graded shape is 
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approximated by multilevel binary steps. Both fabrication processes rely mainly on 
etching techniques that are difficult to accurately control. As a result, the shape and 
depth of the grooves can differ fiom those desired, which leads to reduction of 
diffi*action efficiency and poor repeatability of performance. 
5 In the present invention we consider DOEs based on the space-domain 

Pancharatnam-Berry phase. Unlike diffi-active and refractive elements, the phase is 
not introduced through optical path differences, but results from the geometric phase 
that accompanies space-variant polarisation manipulation. The elements are 
polarization dependent, thereby, enabling multipurpose optical elements that are 

10 suitable for applications such as optical interconnects, beam splitting, and optical 
switching. Recently, we have demonstrated continuous Pancharatnam-berry phase 
optical elements (PBOEs) based on subwavelength grating. However, applying the 
constraint on the continuity of the grating leads to space varying of the local 
periodicity. Therefore, in order not to exceed the Wood anomaly, the elements are 

15 limited in their dimensions as well as optical parameters. Moreover, the results of 
space varying periodicity complicates the optimization of the photolithographic 
process. In order to overcome these limitations, we show that such elements can be 
realized with a multilevel discrete binary geometrical phase, using a computer- 
generated space-variant subwavelength dielectric grating. By multilevel discrete 

20 controlling of the local orientation of the grating, having uniform periodicity, we can 
achieve any desired phase element. We experimentally demonstrate Multilevel 
Pancharatnam-Berry phase Optical Elements (Multilevel-PBOE) as blazed diffraction 
grating and diffractive focusing lens, for the 10.6- //m wavelength from a CO2 laser, 
and show that high diffraction efficiencies can be attained utilizing single binary 

25 computer-generated mask, as well as forming multipurpose polarization dependent 
optical elements. 

The PBOEs are considered as waveplates with constant retardation and space 
varying fast axes, the orientation of which is denoted by Oix^y) . It is convenient to 
form such space varying waveplate using subwavelength grating. When the period of 
30 a subwavelength periodic structure is smaller than the incident wavelength, only zero 
order is a propagating order, and all other orders are evanescent. The subwavelength 
periodic structure behaves as a uniaxial crystal with the optical axes parallel and 
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perpedicutar to the subwavelength grooves. Th«-efore, by fabricating quasi-periodic 
subwavelength structures, for which the period and orientation of the subwavelength 
grooves was space varying, we realized rotating waveplates. Furthermore, we showed 
that such polarization manipulation inevitably led to phase modification of 
5 geometrical origin results firom local polarization manipulation and was in-fact a 
manifestation of the geometrical Pancharatnam-Bcrry phase. 

It is convenient to describe PBOEs using Jones calculus. The space-dependent 
transmission matrix for PBOE is given by applying the optical rotator matrix on the 
Jones matrix of the subwavelength dielectric grating to yield in helical basis, 

where /^^ , are the real amplitude transmission coefficients for light polarized 
perpendicular and parallel to the optical axes, and ^ is the retardation of the grating. 
Thus for an incident plane-wave with arbitrary polarization |E;„ ) , we find that the 
resulting field is, 

15 |E^,> = ^|E,„) + V^e'-<-'lR> + V^e-«*'-'|L) , (2) 

where ;7^=|i(^+f,e")f, t,, =\^(t,-t/*){E^\l.f .»/, =|i(/, -/^'XE. 1^)1'. 

are the polarization order coupling efficiencies , {tx\j3) denotes inner product, and 
I R) = (l O)^ and | L) = (O if represent the right-hand and left-hand circular 

polarization components, respectively. From Eq. (2) it is evident that the emerging 
20 beam fix>m a PBOE comprises three polarization orders. The first maintains the 

original polarization state and phase of the incident beam, the second is right-hand 
circular polarized and has a phase modification of 20{x,y) , and the third is opposite 
to the second at polarization hand and at phase modification. Note that the 
polarization order coupling efficiencies depend on the groove shape and material, as 
25 well as on the polarization state of the incident beam. For the substantial case of 
= = 1 , and ^-n 2Si incident wave with |R) polarization is subject to entire 
polarization state conversion and results in emerging field, 

|E^> = e-"'*''lL>. (3) 
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An important feature of Eq. (3) is the phase factor ^a(^yy)\mod2;r =~2^(^,>')|^2;r ^^at 

depends on the local orientation of the subwavelength grating. This dependence is 
geometrical in nature and originates solely from local changes in the polarization state 
of the emerging beam. 

5 In our approach, the continuous phase function (x, y) is approximated with 

multilevel discrete binary steps leading to formation of PBOE with multilevel discrete 
binary local grating orientation. In the scalar approximation, an incident wave front is 
multiplied by the phase function of the multilevel phase element that is described by, 
exp[/F(^^)], where is the desired phase and F(^rf) is the actual quantized phase. 

10 The division of the desired phase to N equal steps is shown in Fig. 16, where the 
actual quantized phase is given as a function of the desired phase. The Fourier 

expansion of the actual phase front is given by exppFC^^)] = ^jCt ®xp07^rf) , where 
C; is the / th-order coefficient of the Fourier expansion. The diffraction efficiency, rjf , 
of the / th-diffracted order is given by 77, = \Ctf . Consequently, the diffraction 

15 efficiency 7, for the first diffracted order for such an element is related to the number 
of discrete levels Nhy i]^ = {{N 1 7v)s\n{7t I N)f . This equation indicates that for 
2,4,8, and 16 phase quantization levels, the diffraction efficiency will be 40.5%, 
81.1%, 95.0%, and 98.7%, respectively. The creation of a multilevel-PBOE is done by 
multilevel discrete orientation of the local subwavelength grating as illustrated in 

20 Fig.l6. 

Supposing we wish to design a blazed polarization diffraction grating, i.e. a 
grating for which all the diffracted energy is in the 1st order, when the incident beam 
is |R) polarized. We designed a multilevel-PBOE that acts as a diffraction grating by 

requiring that (p^ = {27rld)x\^2n > forming the quantized phase function 
25 F{(pa) depicted in Fig. 16, where d is the period of the diffraction grating. In order to 
illustrate the effectiveness of our approach, we realized multilevel diffraction gratings 
with various number of discrete levels, N = 2,4,8,16,1 26 . The grating was fabricated 
for CO2 laser radiation with a wavelength of IQAfim , with diffraction grating period 
d = 2mm and subwavelength grating period A = 2//m . The magnified geometry of the 
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gratings for nximber of discrete levels, AT = 4 , is presented in Fig. 17, as well as their 
predicted geometrical quantized phase distribution. The elements were fabricated on 
500///W thick GaAs wafers using single binary mask, contact photolithography and 
electron-cyclotron resonance etching with BCI3 to nominal depth of 2.5 fim , resulting 

5 in measured values of retardation (i> = nll, and = = 0.9 . By combining two such 
gratings we obtained a grating with (f> = 7t . These values are close to the theoretical 
predictions achieved using rigorous coupled wave analysis. 

Following the fabrication, the multilevel-PBOEs were illuminated with a 
right- hand circularly polarized beam, [R) , at 10.6//m wavelength. We used the 

10 circular polarized in order to transmit only the desired |L) state, and to eliminate the 
I R) polarization order, which appeared because of the error in the etched depth of the 
grating to achieve the desired retardation of ^ = ;r . Figure 17 shows the measured 
and predicted diffraction efficiency for first diffracted order for the different 
multilevel-PBOEs. There is a good agreement between the experimental results and 

15 the prediction efficiencies, providing verification of the expected multilevel discrete 
binary phase. 

In addition, we formed a multilevel Pancharatnam-Berry focusing element for 
a 10.6// m wavelength, which had a multilevel discrete spherical phase function of 

=F[(2;r/A)(A:^ +/^)''^] , with 1 0 mm diameter, focal length / = 200mm , 
20 and number of discrete levels = 8 . Figure 18 illustrates the magnified geometry of 
a focusing lens based on multilevel-PBOE with = 4 , as well as the predicted 
multilevel geometrical phase. A scanning electron microscope image of a region on 
the subwavelength structure, which we had fabricated, is shown in the inset of Fig. 
16. Diffraction limited focused spot size for |L) transmitted beam was measured, 

25 while illuminating the element with |R) polarization state, and inserting a circular 

polarizer. The inset in Fig. 18 shows the image of the focused spot size as well as the 
measured and theoretically calculated cross section. The measured diffraction 
efficiency was 95% in agreement with the predicted value. The geometrical phase of 
the PBOE is polarization dq^endent, therefore, we experimentally confirmed that our 
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element is a converging lens for incident |R) state, and a diverging lens for incident 
I L) state, as indicated by Eq. (2). Moreover, it is possible to form bifocal lens for 
PBOE with a retardation phase of ^ = ;r while illuminating with a linear polarization 
beam, and inserting refractive lens following the PBOE. Trifocal lens can be created 
5 for PBOE with a retardation phase of ^ = ;t/2 resulting in distinct three different 

focuses for |R), linear, and |L) polarization states. 

To conclude. We have demonstrated the formation of multilevel discrete 
Pancharatnam-Berry phase optical elements using computer-generated space-variant 
subwavelength dielectric grating. We realized blazed diffraction grating, as well as 
1 0 polarization dependent focusing lens. The introduction of space varying geometrical 
phases through multilevel-PBOEs, enables new approaches for novel polarization- 
sensitive optical elements. We are currently investigating a photolithographic process 
with the purpose of achieving accurate control of the retardation phase to yield only 
the desired polarization order. 

15 

5. Near-field Fourier transform polarimetry by use of a quantized space-variant 
subwavelength grating 

Optical polarimetry measurement has been widely used for a large range of 
applications such as ellipsometiy (J. Lee, J. Koh and R. W. Collins, "Multichaimel 

20 Mueller matrix ellipsometer for real-time spectroscopy of anisotropic surfaces and 

films", Opt. Lett. 25, 1573-1575 (2000)), bio-imaging (V. Sankaran, M. J, Everett, D. 
J. Maitland and J.T. Walsh, Jr., "Comparison of polarized light propagation in 
biological tissue and phantoms". Opt. Lett. 24, 1044-1046 (1999)), imaging 
polarimetry (G. P. Nordin, J. T. Meier, P. C. Deguzman and M. W. Jones, 

25 "Micropolarizer array for infrared imaging polarimetry", J. Opt. Soc. Am. A 16, 
1 168-1 174 (1999)) and optical communications (P. C. Chou, J. M. Fini and H. A. 
Haus, "Real-time principle state characterization for use in PMD compensators", 
IEEE Photon. Technol. Lett. 13, 568-570 (2001)). A conunonly used method is the 
measuring of the time-depended signal once the beam is transmitted through a 

30 photoelastic modulator (G. E. Jellison, Jr., "Four channel polarimeter for time- 
resolved ellipsometry". Opt. Lett., 12, 766-768 (1987)) or a rotating quarter wave 
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plate (QWP) followed by an analyzer (E. Collet, Polarized Light (Marcel Dekker, 
New York, 1 993)). The polarization state of the beam can be derived by Fourier 
analysis of the detected signal. This method, however, requires a sequence of 
consecutive measurements, thus making it impractical for real-time polarization 
S measurement in an application such as adaptive polarization-mode dispersion 

compensation in optical commimication. Morieover, it involves either mechanically or 
electronically active elements resulting in a complicated and cumbersome device. 

An increasing demand for faster and simpler methods has led to the development 
of the simultaneously four-channel ellipsometer (R. M. A. Azzam, ^'Integrated 

10 polarimeters based on anisotropic photodetectors". Opt. Lett. 12, 555-557 (1987)). 
Using this method, the beam is split into four channels; each is analyzed using 
different polarization optics while the real-time polarization state is calculated from 
the measured intensities. The main drawback of this method is its high sensitivity to 
statistical errors due to the low number of measurements. Recently, Gori (F. Gori, 

15 "Measuring Stokes parameters by means of a polarization grating". Opt. Lett. 24, 584- 
586 (1999)) proposed real-time polarimetry by use of a space-variant polarizer, in 
what is basically a manifestation of the four-channel technique. This method relies on 
measuring the far-field intensity and therefore is not suitable for on-chip integrated 
applications. 

20 In a recent paper (Z. Bomzon, G. Biener, V. Kleiner and E. Hasman, "Real-time 

analysis of partially polarized light with a space-variant subwavelength dielectric 
grating". Opt. Lett., 27, 188-190 (2002)) we presented a space domain analogy to the 
rotating QWP method. Using this method, spatial intensity distribution analysis is 
applied for real-time, near-field polarimetry. The intensity modulation is achieved by 

25 a space varying wave plate, realized as a computer-generated, space- variant, 

continuous subwavelength grating, followed by an analyzer. However, continuity of 
the subwavelength grating had led to a space dependent local period. Therefore, in 
order for the period not to exceed the Wood anomaly, the elements were restricted in 
their physical dimensions. Moreover, the varying periodicity had complicated the 

30 optimization of the photolithographic process and had led to spatial variations in the 
retardation of the element. 
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In this paper we propose real-time near-field polarimetry by spatial discrete 
rotating of the groove orientation of a subwavelength dielectric grating. The grating of 
this type of element is divided into equal sized zones. The subwavelength grooves are 
of uniform orientation and period at each zone and are rotated at quantized angles, 
5 respective to each zone. The resulting elements are unlimitted in their dimensions and 
have uniform optical parameters. We named such elements quantized space-variant 
subwavelength dielectric gratings (QSGs). Our method is less sensitive to statistical 
errors due to the increased number of measurements, and is suitable for real-time 
applications and can be used in compact configurations. 

1 0 The concept of the near-field polarimetry based on subwavelength gratings is 

presented in Fig. 19. A imifonnly polarized light is incident upon a polarization 
sensitive medium (e.g. biological tissue, optical fiber, wave plate etc.) and then 
transmitted through a QSG, which acts as a space-variant wave plate, followed by a 
polarizer. The resulting intensity distribution is imaged onto a camera and captured 

15 for fiirther analysis. It is shown hereinafter that the emerging intensity distribution is 
uniquely related to the polarization state of the incoming beam. This dependence is 
given by a spatial Fourier series analysis, whereby the resulting Fourier coefficients 
completely determine the polarization state of the incoming beam. 

The QSGs are considered as wave plates with constant retardation and space 

20 varying fast axes, the orientation of which is denoted by Oix.y) . It is convenient to 
form such space varying wave plates by use of a subwavelength grating. When the 
period of a subwavelength periodic structure is smaller than the incident wavelength, 
only the zero order is a propagating order, and all other orders are evanescent. The 
subwavelength periodic structure behaves as a uniaxial crystal with the optical axes 

25 parallel and perpendicular to the subwavelength grooves. Therefore, by fabricating 
quasi-periodic subwavelength structures, for which the orientation of the 
subwavelength grooves is changed along the length of the element, space-variant 
wave plates can be formed. 

The creation of a QSG is done by discrete orientation of the local subwavelength 

30 grating as illustrated in Fig, 20(a). The QSG is obtained by dividing the element into 
N equal sized zones, along the x-axis, where each zone consists of a unifomi 
orientation as well as a uniform subwavelength grating period. The orientation of the 
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15 



grooves is defined by the angle, 9, between the grating vector of the 
subwavelength grating (perpendicular to the grooves) and the x-axis, therefore ^ is a 
function of the x coordinate \6(x) \ , The grating period, rf, is defined as the distance 
between the nearest zones having identical orientations. We consider the period of 
grating d as larger than the incident wavelength X , whereby the local subwavelength 
period of the grooves A , is smaller than the incident wavelength. Fig, 20(a) presents 
a QSG with a period that consists of iV=4 zones of uniform orientation of the 
subwavelength grating. We named N as the number of quantized levels. 

The polarization state within Stokes representation is described by a Stokes vector 
S^(SoSu S2, S3/, where So is the intensity of the beam, whereas 5/, S2. S3 represent 
the polarization state, hi general, 5^ > 5,^ + S2 + S^ , where the equality holds for 
fully polarized beams. The polarization state emerging fix)m an optical system (i.e. 
wave plates, polarizers etc.) is linearly related to the incoming polarization state 
through, S'=MS, where A/ is a 4-by-4 real Mueller matrix of the system and 5, S' are 
the Stokes vectors of the incoming and outgoing polarization states respectively. 

The optical system under consideration consists of a QSG followed by a polarizer. 
This composite element can be described, in Cartesian coordinates, by Mueller matrix 
M=PR('0)WR(0) , (1) 



20 where, Ri0) = 



fl 0 0 d\ 

0 cos 20 sin 20 0 

0 -sin 20 cos 20 0 

0 0 0 1 



is the Mueller matrix that represents 



rotation of the axis fi-ame by angle 0, 



2 



X V 

X y 

0 
0 



0 

0 



0 
0 



0 
0 



2/,/^ cos^ - 2tJ^ sin < 



is the Mueller matrix of a 



2r,r,sin^ 2v,cos^^ 
transversally uniform retarder, wifli retardation ^ and real transmission coefficients 
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10 



for two eigen-polarizations tx, ty. Finally, ^ = 



''I 1 0 

1 1 0 

0 0 0 

0 0 0 



0\ 
0 
0 
0 



is the Mueller matrix 



of an ideal horizontal polarizer. 

The outgoing intensity can be related to the incoming polarization state of the 
beam by calculating the Mueller matrix given by Eq. (1) and using the linear relation 
between die incoming and the outgoing stokes vectors, yielding, 
So'ix) = /.{ASo+ZjiA + QS, +5(5, +S^)cos20ix) + 

(BS^ -DS j)sin20ix)+ (2) 
+ y2{A- C)(S, cos4tf(x) + Si sin4^(x))} , 

where A = tl+tl, B = tl-tl, C = 2tJ^ cos^, D = 2tj^ sin^. 

Equation (2) describes the intensity of the outgoing beam as a truncated Fourier 
series with coefficients that depend on the Stokes parameters of the incident beam. In 
our case 0 represents the quantized rotation angle of the retarder, (subwavelength 
grating), as a function of the location along thex-axis. A single period of dean be 
written explicitly as. 



6'(x) = 



0 

n_ 
N 



N 



N 



0<x< — 
N 

d 2d 
— <x< — 
N 



d(m-\) dm 

— -<x< — 

N N 



n(N-\) d{N-\) 



N 



<x<d 



(3) 



15 where m is an integer number, N^m>0. Figure 2(b) illustrates 0(x) for a number 
of quantized levels A^. Note that 6(xjis periodic in 

Inserting the quantized angle described by Eq. (3), into Eq. (2) and then expanding the 
trigonometric expressions into a Fourier series yields the resulting intensity 
distribution, 
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5o'(*) = X 



, cos- 



iBS2-DS,)f\c„ sin 



2mix . . 2;mxl 



d 

iTtnx 



■ + 6„sin 



2;mx1 
+ cos — — l-h 



r « r 
«=i L 

V. /i-i L 



" d \ 
— +/„sm— J- 

I + «. cos — — 

d " d \ } 



(4) 



(5) 



10 



The corresponding Fourier coefficients in £q. (4) are given by, 

N . (2rc\ 

a. = c. — sm — 

" " 2an \N J 

, , N . 2(^\ 
b„ =d„ = — sin — 

N . f4;r^ ' 

for w = A:A^ ± 1 (4=0,1,2,3,...). and zero otherwise. One can see, from Eq. (5), that 
mcreasing the number of quantized levels, AT, increases a, ci, ei, gi towards unity and 
decreases 6/, rf/,//, hi towards zero. Moreover, as the number of quantized levels 
increases, the higher order terms tend to reach zero as well, thus at the limit of an 
infinite number of quantized levels, Eq. (4) is degenerated to the case of continuous 
space-variant subwavelength grating. 

Using the Fourier analysis, the first coefficients of Eq. (4) yields. 
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(6a) 



4 *2.7DC 




(66) 



(6c) 




0 



(6rf) 



C)(S2^, + S.y; ) = - f So' Wsin(^>fa 



(6e) 



These equations are a linear combination of the Stokes parameters of the incident 
beam. In order to extract S0-S3, Eqs. (6a)-(6e) should represent 4 independent 
equations, which are obtained for > 5 . However, a larger number of quantized 
5 levels, A^, is desirable, in which case a larger portion of the intensity is represented by 
the first harmonics of every series of Eq. (4) and, therefore, a larger signal to noise 
ratio can be obtained. 

Subwavelength quasi-periodic structures are conveniently described as depicted in 
Fig. 20 using Jones calculus. The QSG, which is a birefringent element with optical 

10 axes, (parallel and perpendicular to the grating grooves) that rotate periodically in the 
jc -direction, can be represented as a polarization diffraction grating. When a plane 
wave with a uniform polarization is incident on such a periodic subwavelength 
structure, the transmitted field will be periodic in both the polarization state and 
phase-front, therefore, we can expect this field to yield discrete diffraction orders. The 

1 5 interference between the diffraction orders in the near-field yields spatially intensity 
modulation. Evidently, the resulting interferogram pattern is directly related to the 
incident polarization state. 

In this representation, a uniform periodic subwavelength structure, the grooves of 
which are oriented along the -axis, can be described by the Jones matrix. 



where fj^and ty are the real amplitude transmission coefficients for light polarized 
perpendicular and parallel to the optical axes and ^ is the retardation of the grating. 



20 




(7) 
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Consequently, the space-dependent transmission matrix of the QSG, can be described 
by 

Tcix) = M(^(x))jM-' (0(x)) , (8) 
where ^x) is the quantized local orientation of the optical axis given by Eq. (3) and 
f cosO — sin^^ 

5 Fig. 20(b), and M(0) = \ ^ the two-dimensional rotation matrix. Note 

i^sin^ cos^ J 

that meanwhile, the polarizer is omitted from the optical system, and this will be 
referred further on. 

For convenience, we convert T^(:c) to the helicity basis, therefore, the space-variant 
polarization grating can be described by the matrix, T(jc) = UT^U"* , where 

10 U = "^^^ ' . ! j is a unitary conversion matrix. Explicit calculation of T(x) yields, 

We further adopt the Dirac's bra-ket notations, in which |R) = (1,0)'^ and 
|l) = (0,1)'^ are the two-dimensional unit vectors for right-hand and left-hand 
circularly polarized light. Thus, the resulting field is a product of an incident plane 
1 5 wave with arbitrary polarization |£^„ ) and the space-dependent transmission 
matrix T(x) given by Eq. (9), yields 

\K.)-nA^.)^rj,e'''^^^^^^ , (10) 

where rj^ =i(^^ +//^) , =t(^ -'X^'XEiJL) .^l =i(^ -^^^(E^JR), are 
the complex field coefficients and (^|>9) denotes the inner product. From Eq. (10) it 
20 is evident that the emerging beam from the QSG, which is denoted by | E„„, ) , 

comprises three polarization orders. The first maintains the original polarization state 
and phase of the incident beam, the second is right-hand circular polarized with a 
phase modification of 2d(x) , while the third has a polarization direction and phase 
modification opposite to that of the former polarization order. Note that the phase 
25 modification of the | R) and | L) polarization orders results solely from local changes 
in the polarization state, therefore, is geometrical in nature. 
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Since ^(x) is a periodic fiinction (Eq. (3), Fig 2(b)), the functions e'^^^'^ and 
Eq. (10) can be expanded into the Fourier series. Taking into account the 
connection between the Fourier series of a conjugated functions e'^^*'^ and e"'^^^""^ , 
leads to the equation, 

5 |E^> = ;7.|E.> + '7, |:«.«""""lR> + '7. I:«V^'""'|L> . (11) 

where a„ = r^je'^^^*^e'^'""'^''dx . Based on Eq. (1 1) we find that the diffraction 
efficiency into the m*** difiiacted order of |R) polarization order, ( 7m = )» is 
equal to the diffraction efficiency into the -m* diffracted order of | L) ( 7j!:„, = |cf^ |^ ), 
thus, |R) and |L) diffract adversely. Thereby, as depicted in Fig. 3, once a uniformly 
10 polarized beam is incident upon the QSG, the resulting beam comprises three 

polarization states described by Eq. (10), whereas, the polarization orders of | R) and 

|l) states are split into multiple diffraction orders due to the discontinuity of the 
phase, 2d(x). 

In our near-field polarimetry concept, the wave-front emerging from the 
1 5 polarization grating is incident upon a linear polarizer, described in the helicity basis 

by Jones matrix of the form, P-—\ - The linear polarized field emerging from 
the polarizer is given by the product of Eq. (11) and yielding, 

Consequently, the intensity distribution is given by | , which results after some 

20 algebraic manipulations, to the identical expression as depicted in Eq. (4). 

Note that by using a 4-f telescope configuration for imaging the emerging beam 
from our polarimeter, and inserting a spatial filter in the Fourier plane, the higher 
diffracted order can be eliminated due to the discontinuity of the phase, 20fx). 
Therefore, this can result in an intensity distribution resembling the pattern obtained 

25 by the continuous grating based polarimeter. 



wo 2004/003596 



PCT/IL2003/000541 



44 

The QSG element for CO2 laser radiation of 10.6 ^an wavelength was febricated 
upon 500 fjm -thick GaAs wafer with A - 2//m , d = 2.5mm , and = 16 . The 
dimensions of the element were 30mm x 3mm and consisted of 12 periods of d . 
Firstly, a binary chrome mask of the grating was fabricated using high-resolution laser 
5 lithography. The pattern was then transferred onto the GaAs wafer by use of 

photolithography, after which we etched the grating by electron cyclotron resonance 
with BCU for 35 min, resulting in an approximate 2.5 /^m groove depth. Finally, an 
antireflection coating was applied to the backside of the wafer. Figure 20(c) shows a 
scanning electron microscopy image of a region upon the subwavelength structure, 

10 which we had fabricated, whereas Fig. 2o(d) depicts an image of a cross section of the 
subwavelength grooves. 

Following the measurements, we used the setup depicted in Fig. 19 to test our 
concept for polarization measurements, whereby, a subwavelength metal wire grating 
was used as a polarizer. First, a calibration process was performed by illuminating the 

1 5 QSG followed by a polarizer with three different polarized beams. We determined the 
experimental optical parameters by fitting the curve of Eq. (4) to the measured 
intensity distributions with , , and ^ of the QSG as free parameters. The 

calibration process yielded a QSG having = 0.9, = 0.8 , and retardation 

^ = 0.3;r rad. These values are close to both the theoretical predictions that were 

20 achieved using a rigorous coupled wave analysis (P. Lalanne, and G.M. Morris, 

"Highly improved convergence of the coupled-wave method for TM polarization", J. 
Opt. Soc. Am. A 13, 779-784 (1996)) for a groove profile depicted in Fig, 20(d), and 
to the direct measurement of the optical parameters using ellipsometric techniques (E. 
Collet, Polarized Light (Marcel Dekker, New York, 1993)). 

25 In order to test the ability of our device in conducting polarization 

measurements of fully polarized light, we used a CO2 laser that emitted linearly 
polarized light at a wavelength of 1 0.6 //w and replaced the polarization-sensitive 
medium with a QWP. The images were captured by a Spiricon Pyrocam III, at a rate 
of 24 Hz. Figure 22 shows the measured intensity distributions captured in a single 

30 camera frame, when the fast axis of the QWP was set at angles O*' , 20° , and 45*" , as 
well as the predicted results calculated by Eq. (4) for A=1.45, B=0.17, C=0.8464, and 
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D=l . 165. Consequently, Fig. 23 shows the measured and predicted Stokes parameters 
of a resulting beam as a function of the orientation of the QWP. We determined the 
experimental values of 5, , , and ^3 by using Eq. (6). There is a good agreement 
between the predictions and experimental results. Moreover, Fig. 6 shows the 
5 experimental and theoretical azimuthal angle, {(/ , and ellipticity x > calculated from 
the data in Fig. 5, by use of the relations tan(2^) -S^ISy^ and sin(22') = ^3/^0 [6]. 
The measurements yield a standard deviation error with respect to the theoretical 
prediction of 2.6** and 0.6^ for the azimuthal angle and the ellipticity angle, 
respectively. The errors of the polarization measurements result mainly from 

10 systematic errors such as the nature of the algorithm, imprecision in the rotation of the 
QWP, optical performance of the QWP, low resolution of the IR camera and dynamic 
range, as well as statistical noise due to spatial and temporal fluctuation of the light 
emitted from the laser, shot noise and amplifier noise of the IR camera, pixel response 
nonuniformity, and quantization noise. The standard deviations in the azimuthal angle 

15 and ellipticity for a series of successive measurements were 0.2^ and 0.07** , 
respectively. 

In order to demonstrate the use of a QSG for polarimetry of partially polarized 
beams, two independent CO2 lasers of orthogonal linear polarization states were 
combined by use of the setup depicted in the inset at the top of Fig. 25. The degree of 

20 polarization (DOP) is defined by DOP = ^S^ -¥8^ ^-S^ I Sq . For incoherent beam 

summation, the Stokes vector of the resulting beam is the sum of the Stokes of the 
combined beams. In the case of two orthogonal linear polarized beams, the DOP is 
given by DOP = (/j - l2)K^i + ^2 ) » were I^ and 1 2 are the intensities of the 
horizontally and vertically polarized beams, respectively. Figure 25 shows the 

25 measured and predicted DOP as a function of the intensity ratio /| //2 - The inset 

shows the experimental intensity distributions for two extreme cases. The first is for 
equal intensities (/, = /j), in which the measured DOP is of 0.0059, indicating 
unpolarized light. The second is for illuminating a single laser only (i.e. /j = 0), in 
which the measured DOP results in 0.975, indicating fiilly polarized light. This 

30 experiment shows the ability to obtain all four Stokes parameters simultaneously. 
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thereby, emphasizing the good agreement between prediction and measurement for 
partially polarized light. 

We have theoretically analyzed and experimentally demonstrated the use of a 
computer-generated space-variant quantized subwavelength dielectric grating for real- 
5 time polarization measurement. The quantized subwavelength gratings are unlimited 
in their dimensions, have uniform optical properties (<r, and and in general, are 
lightweight, compact and highly efficient. Both the space-variant subwavelength 
dielectric grating and the polarizer were realized using photolithographic techniques 
commonly used in the production of micro-electric devices. Therefore, the camera and 

1 0 the gratings could be combined into a single chip, resulting in a very small device for 
polarization measurements. 

It should be clear that the description of the embodiments and attached Figures 
set forth in this specification serves only for a better understanding of the invention, 
without limiting its scope as covered by the following Claims or their equivalents. 

15 It should also be clear that a person skilled in the art, after reading the present 

specification could make adjustments or amendments to the attached Figures and 
above described embodiments that would still be covered by the following Claims or 
their equivalents. 



20 



